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Abstract 

A mesoscopic conducting ring pierced by magnetic flux is known to support the persistent electron current. Here we propose 
possibility of the persistent current in the ring made of the perfect crystalline insulator. We consider a ring-shaped lattice 
of one-dimensional "atoms" with a single energy level. We express the Bloch states in the lattice as a linear combination of 
atomic orbitals. The discrete energy level splits into the energy band which serves as a simple model of the valence band. We 
show that the insulating ring (with the valence band fully filled by electrons) supports a nonzero persistent current, because 
each atomic orbital overlaps with its own tail when making one loop around the ring. In the tight-binding limit only the 
neighboring orbitals overlap. In that limit the persistent current at full filling becomes zero which is a standard result. 
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A mesoscopic conducting ring pierced by magnetic 
flux <f> is known to support the equilibrium persistent 
current. [1,2]. Here we propose possibility of the per- 
sistent current in a ring made of the perfect crystalline 
insulator like for instance the intrinsic silicon crystal. 
Consider the ring-shaped lattice of one-dimensional 
(ID) "atoms" in Fig. 1. The electron wave function 
ip(x) in such ring obeys the Schrodinger equation 

h ' " +V(x) 



2m dx 2 



Tp{x) = Eip(x) 



with the cyclic boundary condition [1] 

ip(x + L) — exp(i27r0/ '</>o)ip(x) 



(1) 



(2) 



where m is the free electron mass, x is the position 
along the ring, L is the ring length, (f>o = h/e is the 
flux quantum, and V(x) is the potential of the lattice. 
The persistent current is I — —d E n (cj>) / d(f>, where 



Y2 n En is the ground-state energy of all N e electrons 
in the ring. Solving equation (1) in the lattice model 
with nearest neighbor hopping one gets [1] for odd N e 



I = 



4tt TT sin§N e . 
- U — T J hr- sin 



N<j> 



2?r <j> 



(3) 



where U is the hopping amplitude. The persistent cur- 
rent (3) is nonzero in the conductor (N e < N) but zero 
in the insulator (N e — N). We want to show that the 
persistent current is nonzero also in the insulator. 
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Fig. 1. Potential V(x) in a ring-shaped periodic ID lattice 
of atoms modelled by one-dimensional potential wells of 
width 2d and depth Vb- The lattice period is a. The ring 
circumference is L = Na, where N is the number of atoms. 
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The atomic orbital <p(x) and energy E ai in a single 
isolated atom at x = obey the Schrodinger equation 



h 2 d 1 



if(x) = E at <p(x) 



(4) 



where v(x) is the atomic potential - the same potential 
well as in Fig. 1 but with the infinitely thick barriers. 

We want to express ip(x) vi a the isolated atomic or- 
bitals. We start with a toy model represented by the 
"single-atomic" ring in Fig. 2. In the ring geometry 

V(x + L) = V{x) . (5) 

Due to the condition (5), the ring potential V(x) is 
formally tractable as an infinitely long periodic poten- 
tial with period L (Fig. 2). Therefore, V(x) can be ex- 
pressed as an infinite sum of isolated atomic potentials, 



V(x) 



J2 v ( x - i L ) , 



(6) 



albeit we treat the ring of the finite length L, not the 
infinite crystal. Representation of the infinite crystal 
allows us to expand tp(x) into the atomic orbitals as 



4>(x) 



E 



ikjL 



ip(x - jL) . 



(7) 



Expression (7) can be seen to obey the boundary con- 
dition (2) for k — + n), where n is the integer. 
We set k = ^(4- + n) into the equation (7). We get 



e l2lr *o j <p(x - jL) 



(8) 



j — — oo 

We can now evaluate (i/>|i/>) as 
{ip\tp) — / dxip (x)%j){x) = 

oo oo „ L. 

e 2 ^ (j - j>) / ^ dxtp(x-j'L)<p(x-jL) 

j = — OG j' = -oo ~~ 2" 

°° p OO 

= e' t2w ^^ / dx<fi(x - AjL)(p(x) , (9) 

A . _ J — OO 



where the integral in the second line was rewritten as 
J 2 L jL dxifi(x — [j' — j]L)ifi(x), the variables j and f 

were changed to the variables j and Aj = j' — j, and 
the summation over j was performed. 

Further, we write the Hamiltonian in equation (1) in 
the form H = Hf + V'(x - jL), where 



V(x) 



the ring with 
a single "atom" 
in ID representation 



-3L/2 -L -L/2 



L/2 



3L/2 



Fig. 2. Potential V(x) in the ID ring with a single ID 
atom. The ring circumference L is expanded as a linear 
cell between —L/2 and L/2. The cell is then periodically 
repeated, which is justified by the periodic condition (5). 

is the Hamiltonian in (4) for the atom at site jL and 
V'(x-jL)=V(x)-v(x-jL). (11) 
For the electron energy E — (ip\H\ip) / (ip\ip) we obtain 

1 



OO oo 



E = E at + 



V V e i27r ^ (j - J,) 

j=-oo j' = -oo 



x f 2 dbnp{x - j'L)V'(x - jL)<p{x - jL) , (12) 

where we have applied Hj f ifi{x — jL) = E at ip(x — jL). 
Using similar algebra as in equation (9) we obtain 

2 jr 7A3 cos(27r-f A,) , (13) 



E = E — 7,) - 



<v#> 



where 7Aj = — <p(x — AjL)V' (x)ifi(x) is the over- 
lap integral and we have applied = 7a j • 

We notice that 7a^ <C 7i for |Aj| > 1 and we also 
put (V'l 1 /'} — 1- The formula (13) thus reduces to 



E = E a 



- 70 - 271 cos (27tc?!>/0o) 



(14) 



We set (14) into the expression for the single-electron 
persistent current, / = — dE((f>) /d(j>. We get 



/ = — (47r/0o)7i sin (27r0/0 o ) • 



(15) 



In Fig. 3 we compare the formulae (14) and (15) with 
the results obtained by solving the equation (1) nu- 
merically (the numerical method is described in [3]). 
Indeed, the formulae (14) and (15) work well. The cur- 
rent (15) is nonzero due to the nonzero overlap integral 
71 ; the atomic orbital overlaps with its own tail mak- 
ing one loop around the ring. In terms of hopping, the 
electron hops around the ring back into its starting site. 

In the multi-atomic ring in Fig. 1 the cyclic condition 
(2) holds together with the Bloch condition: 

iP(x+L) = e^^'^ix) , 4>(x+a) = e lka ^{x). (16) 
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Fig. 3. Ground-state energy and persistent current versus 
magnetic flux for the "single-atomic" ring from Fig. 2. Two 
different ring lengths L are considered. The "atomic" pa- 
rameters are Vb = 7.03 eV and 2d = 0.14 nm, they give a 
single atomic level E a t = —3 eV. The formulae (14) and 
(15) (shown in a full line) are compared with the purely 
numerical approach (shown in a dashed line). 

The ring potential in Fig. 1 is formally tractable as an 
infinite lattice potential with period a, given as V(x) = 
Sj=-oo v ( x ~ j a )- This allows to expand ip(x) like 

oc 

$k(x) = J2 e ikja v(x-ja) . (17) 

j=—oo 

Expansion (17) obeys the boundary conditions (16) for 
k — + nj , where n is the integer. 

Proceeding in analogy with equation (9) we obtain 

{lpk n \tpk n ) = I dxtp kn (x)lpk n ( X ) = 
oo oo /" 

y~] y~] e thn{j ~ J )o / dxip(x - j'a)<p(x - ja) 

j — — oo j ' — — oo ~~2 
°° roc 

= N e-' k " A i a dx<p(x - A ja )<p(x) , (18) 

A j = -oo J ~°° 
where the integral in the second line was rewritten as 

— -jo 

J 2 L . dx<p(x — \j' — j]a)ip(x) and we summed over j. 
— — — 

We write the Hamiltonian in equation (1) as H = 
Hf + V'(x- ja), where Hf =-£-J^+ v (x- ja) 
and V'(x — ja) = V(x) — v(x — ja). Using the above 
relations we express E n = {tpk n \H \il>k n ) / (ipk n \ipk n ) as 

E _ E at , 1 V V ik„(j-j')a 

J--OC j'--OC 

x / dx(p(x — j a)V (x — ja)tp(x — ja) , (19) 

T 

where we have applied Hf'i^{x — ja) = E at (p(x — ja). 
Using similar manipulation as in equation (18) and 
approximating (tpk„ \ipk„) — N, we obtain 

oo 

E n = E at -7o-2 J2 7A 3 cos (k n A 3 a) , (20) 

A J= 1 




N = L/a N = L/a 



Fig. 4. Left panel: Persistent current / versus ring length 
L for the ID ring made of the perfect crystalline insulator, 
calculated using the formula (22). THe right panel shows 
the same results multiplied by N 2 to estimate the persistent 
current in a hollow 3D cylinder with cross section of order 
~ L 2 and circumference ~ L. The factor of 1000 appears 
when 10 6 cylinders are measured together. Here we have 
used the parameters Vb = 1.74 eV, 2d = 0.01 nm (giving a 
single atomic level E a t = —0.2 eV) and two different a. 

where 7a 3 - = — <p(x — Ajd)V' (x)tp(x) is the over- 
lap integral and = 7a 3 -- Setting (20) into the 
single-electron current 7 n = —dE n /d(j) we obtain 

/„ = - — ^2 -j±'jA j sm(k n A j a) . (21) 

90 A 3 = l 

We sum I = J2 n In over n = 0, ±1, ±2, • • • ± (N e - 1)/2 
assuming odd N e and spinless system. For N e = N 

I = -(4tt/4> )N 1n sm{2^/4> ) , (22) 

which the persistent current in the insulating ID ring. 
We note that if we keep in (21) only the term with Aj = 
1, the sum I = ~}2 n I n gives the standard formula (3). 

Fig. 4 shows the persistent current (22) as a function 
of L. It decays with L exponentially. Enhancement by 
factor ~ iV 2 can be achieved using the 3D cylinder and 
by factor of 10 3 using 10 6 cylinders on a single chip [2]. 
But even the largest considered L is not achievable 
to present nanotechnology. Due to the ID atoms, our 
model is a minimum model. It shows that the persistent 
current in the insulator exists but it strongly underes- 
timates the effect. For realistic 3D insulators we expect 
larger persistent currents and technologically feasible 
sample dimensions. Another possibility is to fabricate 
the ring-shaped ID Kronig-Penney model at full filling 
which also shows the persistent current [4] . 
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